The INTEGRAL / SPI 511 keV Signal from Hidden Valleys 

in 

Type la and Core Collapse Supernova Explosions 

Hooman Davoudiasl* 1 and Gilad Perez''' 2 

Department of Physics, Brookhaven National Laboratory, Upton, NY 11973, USA 

^Department of Particle Physics & Astrophysics, 

^ ■ Weizmann Institute of Science, Rehovot 76100, Israel 

u 

Abstract 



< 



We examine under what circumstances the INTEGRAL/SPI 511 keV signal can originate from 

decays of MeV-scale composite states produced by: (A) thermonuclear (type la) or (B) core collapse 

supernovae (SNe). The requisite dynamical properties that would account for the observed data 

are quite distinct, for cases (A) and (B). We determine these requirements in simple hidden valley 

models, where the escape fraction problem is naturally addressed, due to the long lifetime of the 

new composite states. A novel feature of scenario (A) is that the dynamics of type la SNe, standard 

' candles for cosmological measurements, might be affected by our mechanism. In case (A), the mass 

CO ■ 

of the state mediating between the hidden sector and the SM e + e could be a few hundred GeV 



(N 



X 



and within the reach of a 500 GeV e + e linear collider. We also note that kinetic mixing of the 



| photon with a light vector state may provide an interesting alternate mediation mechanism in this 



case. Scenarios based on case (B) are challenged by the need for a mechanism to transport some of 
the produced positrons toward the Galactic bulge, due to the inferred distribution of core collapse 
sources. The mass of the mediator in case (B) is typically hundreds of TeV, leading to long- 
lived particles that could, under certain circumstances, include a viable dark matter candidate. 
The appearance of longdived particles in typical models leads to cosmological constraints and we 
address how a consistent cosmic history may be achieved. 
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I. INTRODUCTION 



The recent INTEGRAL/SPI measurement of the 511 keV signal [|IJ from e + e~ annihila- 
tion in the Milky Way reconfirms its puzzling characteristics: it appears that e + e~ pairs 
annihilate at a rate of 1.5 ± 0.1 x lO^s^ 1 in the Galactic bulge, whereas the corresponding 
rate in the disk is 0.3 ±0.2 x 10 43 s . Hence, the bulge-to-disk annihilation ratio B/D ~ 3 — 9 
is implied |Ij]. The size of the signal and the inferred value of B/D make it a challenge to 
explain the data using standard astrophysical sources. 

A leading conventional candidate would be the radioactive /3 + decays of type la Super- 
nova (SN) products, such as 56 Co. These explosions can supply the necessary flux of e + 
and are thought to happen mostly in the central region of the galaxy, with an old stellar 
population [l|]. However, it has been estimated that the fraction of the e + that can escape 
the type la ejecta is not large enough to produce the entire observed signal ||. 

It has also been argued that the bulge rate of SN type la (SNi a ) explosions can be inferred 
by mass-scaling of the rate in early type galaxies that host old stars [|J; this would then lead 
to values of B/D that are too small |IJ (however, this extrapolation ignores the differences 
between the evolutions of the Milky Way bulge and such galaxies). In this case, unless one 
assumes that e + e~ annihilation takes place far (on galactic scales) from the source [|J , this 
ratio cannot be straightforwardly explained. 

Various other sources, such as low mass X-ray binaries ]1|] , have been proposed to explain 
the 511 keV signal. However, generally speaking, all such explanations suffer from rather 
sizable uncertainties and it is not clear if any plausible astrophysical mechanism can account 
for the data. This situation has motivated new physics proposals, such as those based on 
the Galactic dark matter (DM) content; see for example |5| |7|- 

Here, we study a different set of new physics possibilities that do not depend on cosmic 
relic abundances of particle populations, but rely only on supernova production of new 
composite states that decay into e + e~. We consider two distinct scenarios, referred to 
hereafter as cases (A) and (B): 

• In case (A), the new states are produced in type la explosions. 

• In case (B) we consider production in core collapse (cc) supernovae. 

A simple estimate shows why this may, in principle, be energetically possible. Roughly 
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speaking, at a rate of 1-2 explosions per century, O (10 48 ~ 50 ) MeV s -1 of binding energy is 
released into the Galaxy. Obviously, even if a tiny fraction of this energy flux is carried by 
the hidden sector, it would be enough to account for the e + e~ annihilation signal. In fact, 
it was proposed long ago that the 511 keV signal could have been generated by the decays 
of heavy r-neutrinos (now ruled out) produced in cc SN (SN CC ) explosions |Q. 

As mentioned above, the 511 keV data could have reasonable standard explanations ||. 
Here, we do not attempt to refute possible standard astrophysical descriptions of the data. 
However, given the uncertainties involved in these (or other) explanations, we believe alter- 
native proposals merit attention. In particular, the long lifetime of particles in the scenarios 
we consider can naturally address the escape fraction problem present in the conventional 
SN Ia picture. Also, our models could have testable predictions for collider experiments [case 
(A)], and may also provide DM candidates [case(B)], though we do not require that our sce- 
narios give rise to a viable DM particle. We note in advance that, while the new scenarios 
we examine provide advantages over the conventional accounts, they are subject to their 
own challenges and constraints. Nonetheless, they offer concrete and interesting examples 
of how hidden valley dynamics could affect astrophysical phenomena, such as the Galactic 
511 keV signal. 

The 511 keV signal suggests that ~ 10 43 e + e~ pairs are annihilated each second, within 
the bulge at the center of the Galaxy. Astronomical observations require that e + be injected 
into the interstellar medium (ISM) near threshold, at energies below a few MeV per e + [ITU 



We will propose that this injected flux comes from the decay of a new state X of mass 
mx ~ 1 MeV, from a hidden dynamical sector of the type discussed in Ref. [11]. We 
further assume that X only decays through SM final states and its branching ratio into 
e + e~ is typically 0(1). Given the complexity of the astrophysical system and the fact that 
the hidden dynamics is non-perturbative our analysis is only done in a semi-quantitative 
manner. Yet, the analysis presented below is concrete enough to make the point that 
hidden valley models may play an important role in dynamics of astrophysical systems and 
the cosmological evolution of the universe. 

In case (A), we find that a successful scenario based on SNj a requires a mediator of mass 
0(300) GeV, connecting the SM and the hidden sectors, making its discovery at a 500 GeV 
e + e~ collider a possibility. We also show how present experimental constraints may be 
addressed. Here, we also briefly discuss how kinetic mixing of the photon with a light vector 
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can provide another potential mediation mechanism between the hidden and the SM sectors 
fT2| . The production of the X particle is Boltzmann suppressed, hence it is not relativistic 
and its e + e~ decay products have MeV-scale energies, as required by observations. The 
lifetime of X is long enough to exit the SN ejecta before decaying and avoids the conventional 
SNi a positron escape fraction problem. The SNi a Galactic population is centrally peaked 
and helps explain the spatial distribution of the 511 keV signal. Potentially troublesome 
long-lived particles produced in the early universe can be efficiently annihilated via strong 
dynamical processes within the hidden sector. A major concern in this scenario, apart 
from direct constraints on the model, is whether Big Bang Nucleosynthesis (BBN) can go 
through without much perturbation. Although we show that this may be plausibly achieved 
within the considered framework, perhaps a more detailed analysis is required to reach a 
firm conclusion. 

In case (B), the X production in SN cc is not Boltzmann suppressed and can be brought 
down to the required levels by lowering the coupling of the hidden and SM sectors. This is 
achieved by assuming a high SM-hidden mediation scale of (9(100) TeV, which leads to a 
long lifetime for X. This feature allows the e + e~ decay products to be deposited far from the 
SN CC sources that are mostly found in the Galactic disk. This long decay length could lead to 
acceptable B/D values, since it causes the surface brightness of the disk e + e~ annihilations 
to be reduced and may also help with magnetic transport of some positrons into the bulge. 
Even though the kinetic energy of the initially produced hidden "partons" is O(10) MeV, 
fragmentation into X and other hidden hadrons softens the final state e + e~ energy to the 
desired few MeV level. Here, acceptable cosmic abundance of long-lived particles may require 
low reheat temperatures, but simple models give rise a to a plausible DM candidate. The 
main shortcoming in this scenario stems from the fact that the disk population of the SN CC is 
bigger than that of the bulge and a Galactic mechanism is needed to transport the generated 
positrons into the bulge, as we will discuss. 

In the following section, we discuss the microscopic structure of our models. In sec- 



tion |rn| , we explain how the 511 keV signal can be addressed via type la supernova explo- 
sions. Section [TV] deals with the core collapse case. In section [V], we discuss the cosmological 
implications of our setup. Section |VT| contains our final conclusions. 
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II. MICROSCOPIC SETUP 



We assume that the hidden sector has a dynamical scale ~ 1 MeV at which the new 
states emerge. This scale sets the masses of the new hidden hadrons 1 . The decays of the 
composite states then yield the requisite e + input at C(MeV) energies . Of course, if X 
is emitted at energies much larger than its mass, the positrons will be very relativistic, in 
conflict with observational evidence that demands a cold positron flux. However, below we 
will argue that Boltzmann suppression, in case (A), or parton shower and fragmentation, 
in case (B), makes this unlikely, given the relevant energies. We will denote the fermionic 
"quark" degrees of freedom in the hidden sector, subject to non-trivial hidden dynamics, by 
/. Below the scale A^, strong dynamics leads to the emergence of hidden hadrons made up 
of / quarks. 



Following Ref. |TT], we will assume that the hidden dynamics is coupled to the SM via 
higher dimension operators suppressed by a scale M. For concreteness, we will consider the 
following dimension-6 operator that can arise in a variety of models: 

(jj i T 1 ij j )(f a T 2 f b ) 

M 2 ' 1 } 

where ipi is an SM fermion and r 12 represent a product of Dirac matrices. For example, if 
the visible and the hidden sectors couple through a massive scalar or vector $ of mass m$, 
then one gets M ~ m$ / ^g\g 2 Pi , where g\ and g 2 denote couplings of the $ to the SM 
and hidden sectors, respectively. In what follows, we focus on the case where $ is a vector 
particle. In our discussion of case (A), we will briefly consider how the kinetic mixing of a 
light vector particle with the photon [0] can provide another mediation mechanism between 
the hidden and the SM sectors. 

In places where a concrete model is appropriate to consider we will adopt the One-Light- 



Flavor (1LF) model discussed in Ref. JTT|. We will denote the lightest pseudo-scalar by r/h 



and the lightest vector by tOh, with m Uh ~ JnJ m Vh based on the results from Ref. fll3| ; 



is the number of hidden colors. There is also a scalar, a^, that is typically intermediate in 



mass and lifetime |L1] . The generic expectation is that cr^ is up to 30% heavier than r)h [11 



1 For simplicity, we do not consider the possibility of having very light pseudo-Goldstone bosons (PGBs) 
well below A^. Inclusion of PGB states would further enrich the structure of the theory and would add 
new scales to the system, hence it is worth future study. 
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and hence it would not decay into r]h + e + e~ , within our typical parameter ranges. Then 
the key aspects of the phenomenology we are interested in can be captured by considering 
only rjh and Uh, as we do for simplicity in the following. We also briefly discuss below the 
possible implications of the — > r/h77 loop-mediated decay. 



III. CASE (A): TYPE IA SUPERNOVA EXPLOSIONS 

Let us begin with an estimate of the X-production in SN Ia . These events are attributed 
to the explosion of accreting white dwarfs through thermonuclear reactions that disrupt the 
star entirely, releasing ~lx M « 2 x 10 33 g of binding energy in the process. Observations 
of similar processes in other galaxies, including those of novas that are believed to share the 
same progenitor, point to a centrally peaked distribution of SNi a It is then plausible to 
assume that they could give rise to B/D > 1 , as long as one can provide a way for a large 
enough fraction of the generated e + flux to escape the explosion and reach the interstellar 
medium (ISM). We will base our estimates on the carbon deflagration model of SNi a , a 



detailed discussion of which can be found in Ref. [fH] ; we will adopt the parameters of the 
widely used W7 model. 

The SNi a explosion occurs over about 1 s, during which a large fraction of the white 
dwarf mass, roughly 0.7 M & , reaches temperatures of order Ti a « 6 x 10 9 K. We will assume 



a mean density p ~ 3 x 10 8 gcm 3 [|I4| . For the electrons in the white dwarf, bremsstrahlung 



processes in which an electron scatters from the electromagnetic filed of the nucleus N are 



the dominant process: eN — > eNX [15|. The rate for scalar or vector emission from non- 



relativistic and non-degenerate plasmas has been calculated in Ref. fL6|. Since Ti a < m e , and 



degeneracy is mild for our reference parameters |T5"|, we adapt their result for a reasonable 



estimate of energy loss rate, in erg g 1 s 1 , |T3|, [L6| 

s x ^e- mx / T ^a'r l 2.8 x 10 26 T 8 a5 F e p]T (2) 

where a' = g /2 /(47r), with g' the coupling of the scalar (massive vector) X state to the 
electron, corresponding to rj = 1(3). Here, Tg = T/(10 8 K), Y e ~ 0.5 is the electron number 
fraction relative to baryons, Xj is the mass fraction of species j with atomic number Zj and 
mass Aj] we will take ^ZjXjZ 2 - j Aj ~ 14 for our estimates. In the above formula, we have 
introduced a factor for Boltzmann suppression, since Ti a < mx- We require that the above 



rate yield ~ 10 43 X particles per second, in order to generate the inferred e + e~ annihilation 
signal. Taking the SNi a rate to be roughly 1 per century in the Galaxy and mx = 2.0 MeV, 
we get 

«'~4.6xl0" 22 /?/. (3) 

With i] = 3 for a vector X and a naive estimate g' ~ m x /M 2 , we find M ~ 300 GeV; 
here and throughout this work we will assume that the decay constant of a hidden hadron 
fx ~ rn x - 

A rough estimate of the X lifetime in this scenario is 



mxg' 2 m\ vauu^eV/ V mx 



16tt 16ttM 4 _ / M \ 4 /2MeV\ 5 



The size of Tx also provides a way to address the problem of e + escape fraction. To see this, 
we argue that the X particles will typically fly out with a velocity vx, much faster than 



the stellar explosion which moves at a relatively low speed t>i a ~ 0.03 \Ti\. The reason for 
this is that the energy distribution of the emitted X particles inherits the electron energy 
distributions inside the SN Ia , hence we expect 

«x~(T, a /m I ) 1 / 2 » Wla . (5) 

Finally, we note that the size of the progenitor of a SNi a is given by the radius of the 
accreting white dwarf, R w ^ ~ 10 4 km ~ 0.03 s. Then, Eq. (Q) implies that the X particles 
will decay well outside the stellar explosion front and will not be absorbed by the ejecta. 
This avoids the escape fraction problem of the conventional SNi a picture, where the positrons 
can get absorbed efficiently as they are released within the outward moving stellar matter. 



A. Direct constraints 



The mass scale M may seem unacceptably low, in light of the existing LEP data |17-19| , 



or direct detection data pOfl . However, we argue that simple assumptions about the nature 
of the coupling of the hidden and SM sectors to the mediating particle $ (such as a heavy 
Z') can allow for such values of M. In particular, if we assume that the hidden sector 
couples to $ strongly with gi ~ 10, whereas the SM couples to $ weakly with g\ ~ l/g2, 
then many existing bounds can be satisfied for m$ ~ 300 GeV. As we do not require X 
production from SM hadrons in SN Ia , a simple assumption would be to take the Z' coupling 



to quarks to be negligible, in which case the Tevatron bounds [2(| would not apply. Here 
we would like to mention that kinetic mixing of the photon with a light vector particle can 
offer another interesting possibility. For example, if the mixing is governed by a loop-level 
mixing parameter e ~ 10~ 3 , for g 2 ~ 0.1 and m$ ~ 1 GeV, we may achieve the required 
level of coupling between the hidden and the SM sectors. Such parameters can accommodate 
precision measurements of (g — 2) e pif , while potentially accounting for a possible deviation 
of (g — 2) M from the SM prediction |p2| . However, for concreteness we will concentrate on 
the case with m$ ~ 300 GeV, in what follows. 

The LEP bounds on (e + e~) v+A — > (e + e~) v+A , from Ref. fl7fl , roughly require A > 7 — 
9 TeV (depending on the sign of the effective dimension-6 operator), where A = \/2tt m$/ g\. 
For gi pa 1/10 and m$ 300 GeV we then get A « 8 TeV which is in the realistic regime. 
We hence adopt a minimal case where $ couples dominantly to (e + e~) v+A in the SM. We 
note that this would introduce a suppression by a factor of 1/2 in Eq. (|2]), due to the random 
polarization of the electrons in the star. Even then, M ~ 300 GeV still yields about 7 x 10 42 
e + e~ pairs per second which, within the accuracy of our analysis, is of the right size. 



Another possible LEP bound comes from e + e — > 7 + fi [|T^, 19" , where the missing 
energy is carried by the invisible // jets. We estimate the cross section for e + e~ — > 7// by 

I a n h c E 2 cm 
e/7 ~ 4 X 7T X 12ttM 4 ' 



(6) 



where the factor of (1/4) accounts for chiral couplings of the initial and final states, the 
factor of a/n accounts for the initial state radiation, and E cm is the center of mass energy. 
In deriving the above formula, we have simply adapted similar results, say, for e + e _ — > 
[HI- A n analogous process arises in models with large extra dimensions 0], where 



the missing energy is carried off by Kaluza-Klein gravitons p5| , |26| . Using the results of 
Ref. [26], we then infer that M 300 GeV is a safe suppression scale, for E cm pa 200 GeV 

on e + e~ — > 7 + fj . 



and n h c = 2, given the LEP2 bounds [18, 19 



B. Astrophysical constraints 

The above inferred value of a' is not constrained by astrophysical processes at T m^, 
such as those of red giant cores at T ~ lOkeV. However, SN CC explosions, characterized 
by T ~ 30MeV, could overproduce the 511 keV signal and, in principle, lead to severe 
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constraints on our SNi a -based scenario. Remarkably, for M ~ 300 GeV, we roughly infer 
a mean free path A ~ M 4 /T 5 ~ 10 2 m for the hidden sector / quarks (even though we 
are assuming that the hidden sector has negligible couplings to SM hadrons of the core, we 
may expect that such a hot environment provides a photon and e + e~ thermal bath that 
can efficiently generate a population of /-quarks). Since the size of the hot SN core is 
10 — 100 km, we see that the / quarks are quite likely trapped and can only be surface- 
emitted. This avoids constraints coming from over-cooling of the SN CC core, which would 
otherwise halt the explosion. 

One may also worry that a trapped "hidden-sphere" will contain a substantial fraction 
of the explosion energy and overproduce the 511 keV signal. However, given the estimation 
of the X lifetime in Eq. (|]) we claim that over-production is not a worry. The reason is 
that the progenitors of the SN CC are red or blue giant stars of radius ~ 10 2 R Q ~ 10 2 s, much 
larger than the decay length of X. Hence, the e + e~ decay products of the X particles will 
be absorbed well- within the SN CC progenitor and will not contribute to an over-population 
of e + in the ISM. In fact, this effect might help to prevent the shock wave (which drives the 
cc explosion) from being stalled, a major outstanding problem in understanding SN CC . 

Regarding cx^, we note that it is expected to be at most 30% heavier than r\ h Hence, 
within our typical parameter range, it would decay into rjh 77, via a one-loop radiative 
process, with an approximate life time 



where Am/, = m ah —m Vh ~ 0.4 MeV. This is roughly consistent with the bound Jl5| produced 
by the Solar Maximum Mission (SMM) satellite |27] looking for a flash of photons emitted 
from the 1987a SN CC , yet may require Amj to be a factor of few smaller. 

IV. CASE (B): CORE COLLAPSE SUPERNOVAE 

As noted before, SN CC can also produce the requisite e + flux to account for the size of the 
511 keV signal. However, as we will see, this can be done in a very different regime of hidden 
sector models, due to the much higher temperature of such explosions. The estimated rate of 
energy release from the SN CC , for particles that free stream out of the core, can be obtained 
from that of vv emission from a non-degenerate stellar core, via NN — > NNuu, with N a 
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nucleoli [H5 



e„ = 2.4 x 10 17 erg g" 1 s" 1 p 15 T* e 5 v . (8) 
Here, the density p 15 is in units of 10 15 g/cm 3 and the temperature T Me v is i n MeV. In 



the SN cc core, T ~ 30 MeV and the rate saturates at pis ~ 0.03 fL5J , and we find e v ~ 
9.6 x 10 24 ergg -1 s -1 . The temperature of the core drops over a time scale of order a few 
seconds. Given the steep temperature dependence of the rate in Eq. (H), we estimate the total 
/-quark (hidden sector) output from the emission over the first few seconds, by rescaling 
the Fermi constant Gf — 1-2 x 10~ 5 GeV~ 2 of the SM weak interactions. For an object 
of mass ~ 1.5M ~ 3 x 10 33 g, the ratio Rf of /-quark output to the total energy budget 
~ 3 x 10 53 erg is given by 

R f ~ 10 4 (M 2 G F )- 2 . (9) 

Based on our preceding discussion, Rf ~ 10 -6 is a reasonable estimate for the required 
output. We then get M ~ 100 TeV as a natural scale for the operator in Eq. ([I]), if one 
wants to explain the 511 keV signal from SN CC production. 

Once the /-quarks are produced, they would "hadronize" on their way out of the SN 
core. Showering and fragmentation of the /-quarks will then distribute their emission energy 
among several X and other hidden hadrons. As a rough guide, we will consider the average 
multiplicity of QCD hadrons in e + e~ collisions. Given the typical SN temperature of 
about 30 MeV and A^ ~ 1 MeV, the data at ^/s ~ 30 GeV will provide a good analogue, 



and we find Nh ~ 20 |p8fl . We then expect O(10) hidden hadrons to result from /-quark 
emissions at T ~ 30 MeV. As long as the initial emitted parton fragments into at least a few 
leading hadrons, their energies will typically be below ~ 10 MeV, which will eventually yield 
positrons at few MeV energies, as required by the data. This feature is a universal aspect 
of generic strong dynamics, regardless of the specific model. Looking at a compilation of 
measurements for center of mass energies between 14-44 GeV, by the Tasso collaboration Jj9| . 



we find, indeed, that in e + e~ — > qq processes less than 1% of the events contain a particle 
which carries 50% or more of the beam energy and the majority of final state hadrons are 
soft. 

A reasonable fraction of these hidden sector hadrons must decay well-outside the 
SN CC progenitor star, if one wishes to disperse the e + e~ throughout the ISM. Since the 
typical radii of the collapsing stars is O(10 2 R Q ), we require Tx ~ 100 s. Note that the naive 
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expectation for the lifetime of X in this scenario is now 

16ttM 4 inll ( M \ 4 /2MeV\ 5 

which is more than sufficient to satisfy this requirement. As we will discuss below, this long 
lifetime could be a useful feature in explaining the distribution of the 511 keV signal. 

The inferred e + e~ flux seems to come mainly from the Galactic bulge and is not evenly 
distributed throughout the Galactic disk. However, the SN CC events are expected to be 
mostly concentrated in the disk. The Galaxy can be roughly described as a disk of radius 
tq ~ 15kpc. We approximate the central bulge of the Galaxy as a sphere with a radius 
tb ~ 3kpc. A reasonable lower bound on the fraction of the cc events that come from 



the bulge region is obtained by the ratio of the disk area inside and outside the bulge |30| 

fB>r 2 B /(r 2 G -r 2 B ), (11) 

which yields fs ~ 0.04. One can also get estimates of the fs using observational data on 
neutron stars and pulsars [[H] . The distribution n cc of SN CC in our galaxy can be described, 
in galactocentric cylindrical coordinates (r,z,9), by pH], |32| 



n cc (r) oc r^e r ^ u 



0J9 e -(,/0.212)" + _ 21 e -(*/0.636)= 



(12) 



where £ and u are fit parameters, and z is in kpc. The neutron star data yield £ = 4 and 
u = 1.25 kpc, whereas the pulsar distribution implies £ = 2.35 and u = 1.53 kpc |3~Ifl . We 



find that the first set of parameters yields fs ~ 0.04, whereas the second set gives fs ~ 0.11. 
In either case, we need to address the paucity of bulge sources. 

Ref. [|[| has suggested that if positrons can propagate away from the Galactic disk they will 
fill out a larger volume and reduce the surface brightness of the disk annihilation emission. 
Also, once the positrons get out to z > 3 kpc, they will escape beyond the cosmic ray 
halo (CRH), where Galactic (poloidal) magnetic fields may transport them into the bulge. 
Ref. argues that, if these conditions are satisfied, then B/D as low as ~ 0.5 can be 
consistent with the current INTEGRAL/SPI observations. However, a naive estimate of 
the escape fraction / esc from the CRH, using slow-down and confinement time scales for 
positrons, yields / csc ~ 0.1 [|J. This is not sufficient to get the required level of disk e + 
population into the bulge, yet there may be room for enhancing / esc , given the uncertainties 
involved. We will next discuss how in our scenario / esc 3> 0.1 can be naturally obtained. 
This would easily allow us to take advantage of the mechanism suggested in Ref. ||. 
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The long lifetime Tx in Eq. (|T0|) is helpful in releasing the requisite e + e~ away from 
the cc explosion and the Galactic disk. Here, we give a rough estimate for the fraction of 
the e + e~ trapped within the CRH. We only consider e + e~ from decays of the X particles 
in our scenario. For simplicity, we approximate the CRH by a slab of thickness d and 
assume that the planar extent of the slab (Galactic disk) is large compared to other length 
scales. Furthermore, it is also assumed that the SN CC explosions are confined to a thin layer 
(set by the thickness of the Galactic disk) in the middle of the slab, which is a reasonable 
approximation. The fraction of the X particles that decay within the slab is at most of 
order 

F x ~- r /2 d6(l-e- 1 ^), (13) 



7T JO 

U2 , 



where I = d/ (2 cos 6). Guided by Eq. fliCf), let us assume that t x ~ 10 s, equivalent to a 
decay length of order lOkpc; this corresponds to nix ~ 2m e . For d = 6kpc, we then find 
F x ~ 0.5; for r x equivalent to 20kpc (M slightly larger than 100 TeV) we get F x ~ 0.3. 
Hence, we see that for typical values of the lifetime Tx, we can easily suppress the fraction 
of the X particles that decay within the CRH to F x ^5 0.5 . The long distance decays of X, 
compared to the Galactic disk scale, ensures that the e + e~ is dispersed over a larger volume, 
which helps diminish the surface brightness of the disk e + e~ annihilation emission [|J. In 
addition, this provides a natural mechanism for transporting the e + e~ beyond the CRH, 
where they could be guided into the bulge by the large scale poloidal magnetic field. Hence, 
the success of case (B) in explaining the 511 keV signal depends on whether enough positrons 
can be transported into the bulge. Whether or not this condition is satisfied in our Galaxy 
is a question beyond the scope of our paper and poses a challenge to a scenario based on 
SN CC . 

The above discussion outlines how hidden dynamics of the type discussed in Ref. JTT 



can in principle produce the gross features of the observed 511 keV signal. However, some 
questions, such as that of the consistency of this type of scenario with standard cosmology, 
can only be answered within more specific models. As a simple example, we will adopt the 
One-Light-Flavor (1LF) model discussed in Ref. [|ll] and mentioned above in section |H[ The 
key aspects of the phenomenology we are interested in can be captured by considering only 
rjh and Uh and not o^, as we do for simplicity in the following. 
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V. COSMOLOGY 



In the early universe, once the cosmic temperature T falls below the dynamical scale 
A/j ~ 1 MeV, various hidden hadronic states appear. The fast interactions amongst the 
hidden sector hadrons allow them to annihilate into the lightest state rjh, as the universe 



gets cooler [11]. If the lightest state has a sufficiently short lifetime, r < Is, then the 
hidden sector decays into the SM degrees of freedom in time for BBN. However, this is 
generally not the case in our type of scenarios. Also, a post-BBN hidden hadronic gas could 
be unacceptable, as such an ensemble reshifts like matter and would come to dominate the 
universe well-ahead of the standard epoch near T ~ 3 eV. Let us examine the status of the 
1LF model regarding these questions. From the results in Ref. []TTj| , the lifetimes of rjh and 
Uh (assuming that the same type of expressions are valid for m ~ 1 MeV) have the following 
dependencies 

and TuhKM*/™^, (14) 
respectively; as mentioned before we will take /„ ~ m„. . 



A. Case (A) 

With M « 300 GeV, n h c = 2, and m Vh;U j h ~ 1.5,2 MeV, respectively, we find r Vh ~ 10 21 s 
and T UJh ~ 3 s. Here, Uh has the kind of lifetime needed for X in case (A). However, r]h is not 
sufficiently long-lived to be a safe DM candidate |33) and if produced in equilibrium, could 
come to dominate over radiation during the BBN, which is not an acceptable outcome. To 
examine these questions, let us estimate the decoupling temperature of the hidden sector. 

Above T ~ mx, we can treat / quarks as free, and their production in the plasma is 
governed by their interactions with e + e~. The rate for this interaction is given by r e f ~ 
n e a e {V. Here, n e ~ T 3 , the cross section a e { is given by Eq. (||), but without the factor of 
a/ir, the relative velocity v ~ 1, and E CTa ~ 2T. The process e + e~ — > ff decouples when 
r cf < H, with the Hubble rate H ~ 1.7 gl /2 T 2 /M Pl ; g* is the number of relativistic degrees 
of freedom and Mpi ~ 1.2 x 10 19 GeV is the Planck mass. The decoupling temperature, for 

~ 10.8 (T > mx) is then given by w 4 MeV. Here, is in principle high enough 
for successful BBN. A similar calculation also suggests that the X particle, identified in 
our 1LF example as Uh, would be out of equilibrium with the SM plasma for T < mx- A 
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simple assumption would be that the reheat temperature was in the 1-4 MeV range, such 
that the /-quarks did not equilibrate with the SM 2 . However, arranging for such a low reheat 
temperature would require a rather non-standard cosmic evolution (for bounds on the scale 
of low temperature inflation see e.g. and references therein). We will next argue that 
the strong interactions amongst the hidden sector states could efficiently suppress the relic 
density of rjh through the relatively fast decays of Uh, in equilibrium. 

The r)h — u)h system stays in equilibrium through fast hadronic interactions governed by 
a cross section ov ~ l/m 2 ;j . As the Uh population decays, its number density decreases as 

n(T) ~ T 3 e~ r/H , (15) 

where r ~ 2.1 x 10~ 25 GeV is the width of ujh, for m u]h ~ 2 MeV. The fast thermal interactions 
decouple once nav ~ H, which yields ~ 0.1 MeV, and hence 

n{T d )/Tl ~ e - r / H ^ ~ 1(T 20 . (16) 

This is a very rough estimate, but shows that the strong hidden dynamics can efficiently 
suppress the rjh number density to negligible levels. We note that a more detailed analysis 
may be called for to determine whether the proximity of the BBN era and the onset of 
hadronization and decays in the hidden sector does not cause large deviations from the 
standard Big Bang picture. It is important to note, however, that the X particle density just 
below the MeV temperature, when BBN starts, drops very rapidly with falling temperature, 
since n(T) oc e~ r ^ H and Y/H oc m x /T 2 . For instance, for rrix = 3 MeV (still in the right 
range to produce enough X particles from SNj a ) and T = 1 MeV we find that n is already 
a few percents of its original density which, in this small perturbation to the 

cosmological energy density. 

B. Case (B) 

For M w 100 TeV and m VhjUh ~ 1.5, 2 MeV, we get T Vh ~ 10 31 s and ~ 10 n s. Thus, 
given our preceding discussion, we see that Uh can easily be a good candidate for particle 
X, responsible for the SN CC -generated 511-keV signal. Note also that T Vh is consistent with 

2 A low temperature mechanism for baryogenesis |34], [55| could provide the baryon asymmetry of the 
universe. 
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511-keV flux bounds [33fl, even if t]h is as abundant as DM. However, in this case rjh will 
be very long-lived and could upset the standard picture of cosmology if its relic density is 
too large. In case (B), the hidden sector decouples from the SM in the early universe at a 
temperature « 13 GeV, where we have assumed ~ 86. Again, the simplest assumption 
is that inflation resulted in a low re-heat temperature of order (within the SM sector) 
that did not lead to any significant production of long-lived hidden hadrons. Nonetheless, 
given the longevity of rjh, it is interesting to see whether there is a reasonable scenario in 
which this particle is the dominant DM of the universe. We will examine this possibility 
next. 

Let us assume that the primordial SM and hidden sectors evolve in a decoupled fashion, 
as discussed above, such that they do not come into thermal equilibrium. In this case, once 
the hidden sector's temperature Th falls below T* ~ 1 MeV, its constituents hadronize. All 
the heavier hidden hadrons will quickly annihilate down to rjh which will then redshift as 
matter. Let us denote the temperature of the SM sector at this point by Tj. In order for 
matter-radiation equality to take place near its standard temperature of Tf ~ 3 eV, we then 
demand 

T*(Tf/Ti) 3 ~ Tj , (17) 

which yields 

10 MeV < Ti < 10 2 MeV for 0.3 MeV < T* < 1 MeV. (18) 

We note that the hidden dynamical scale corresponding to T* can be somewhat lower than 
the hadron masses (as in QCD). Also, over the above range of T i; the SM and hidden sectors 
stay decoupled, as Td ~ 13 GeV. Hence, in this setup rjh can be a realistic DM candidate, as 
long as the conditions ( |18D are met. To complete this discussion, we then suggest a scenario 
in which these conditions can be realized. Imagine that the universe went through a period 
of high scale inflation with a reheat temperature T r 3> Tj, giving rise to a thermalized 
plasma of SM and hidden species. However, a late and milder inflation can cool this plasma 
down to T ~ 1 MeV. As long as the late inflaton only decays into SM degrees of freedom 
and leads to a reheated SM plasma at T = Tj, consistent with flTSD , rjh will survive as a 
viable DM particle. 



15 



VI. DISCUSSION AND SUMMARY 



We examined the conditions under which the 511 keV signal could be related to supernova 
production of an MeV-scale composite state X of hidden dynamics in simple models. Our 
proposed setups naturally avoid the escape fraction problem of conventional SN Ia scenarios. 
However, the models we considered face some challenges, from cosmology or in relation to 
the spatial distribution of the signal. Nonetheless, these proposals offer interesting examples 
of how hidden valley dynamics could affect astrophysical observations, such as that of the 
511 keV signal. Our analysis is done at a semi-quantitative level and a much more detailed 
study is required in order to make definite numerical statements. 

Nevertheless, our study demonstrates, in a concrete manner, that hidden valley dynamics 
could play an important role in cosmology and astrophysics, while also having implications 
for discoveries at accelerators. We considered how SNj a [case (A)] and SN CC [case (B)] explo- 
sions can generate new composite states with relevance to the 511 keV data. 

Case (A): The SN Ia explosions tend to be concentrated mostly near the central part of the 
Galaxy, where the signal originates. Conventional explanations based on SN Ia suffer from 
the small escape fraction of e + generated in radioactive decays of the ejecta. However, in our 
scenario, the X particles have a large enough lifetime to release the e + e~ decay products 
away from the relatively slow moving explosion front. Boltzmann suppressed production 
of X ensures a non-relativistic emission process and leads to a soft e + flux. A minimal 
setup requires a state <3> mediating between the hidden and SM sectors, with a mass m$ ~ 
300 GeV, weak coupling to e + e~ and strong coupling to hidden /-quarks. In this case, $ 
may be directly discovered at a high luminosity e + e~ TeV-scale collider. In addition, in this 
case, it is interesting that both type la and core collapse supernova dynamics are modified 
which may lead to other future signals in precision measurements. If the above effect has 
some z dependence it might lead to systematic effects for the next generation supernova 
observatories. Economical hidden valley models can capture the generic requirements to 
produce the 511 keV signal and lead to acceptable cosmic scenarios, though a more detailed 
analysis may be required to examine the latter. 

Case (B): An alternative possibility is offered by SN CC . These explosions are characterized 
by temperatures far above mx and can easily provide the required X flux. The X production 
is through hidden sector quark jet fragmentation. This results in a soft spectrum for the 
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final state e + e~, as required. The requisite hidden sector production is governed by a large 
suppression scale m$ ~ 100 TeV, well beyond foreseeable collider reaches, which leads to 
long-lived states. In this case, simple hidden valley models may potentially include a DM 
candidate. Given that most SN CC explosions take place in the Galactic disk outside the bulge, 
the e + e~ produced in this case must be transported away from the source. Achieving the 
requisite transfer of e + away from the disk over such length scales may be a challenge within 
a conventional setup. However, in our scenario, the implied decay length of X is typically 
several kpc, allowing for an efficient deposition of e + e~ far from the disk. This leads to a 
reduced disk surface brightness by diffusing the e + e~ over a larger volume and may help 
transport of e + into the bulge, via large scale magnetic fields of the Galaxy. These effects 
could help explain the B/D ratio inferred from observations. However, the relevance of case 
(B) to the 511 keV signal depends on whether the aforementioned transport mechanism is 
operative in the Galaxy. 
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